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We present observation of large, expanding string loops forming around a heated wire tip embed- 
ded in a nematic liquid crystal sample. Loops expand due to convective stretching. This observation 
leads to a new insight into phenomena which could occur in the early universe. We show that local 
heating of plasma in the early universe by evaporating primordial black holes can lead to formation 
\^ • of large, expanding cosmic string loops, just as observed in the liquid crystal experiment. Intercom- 

' mutation of string loops from neighboring black holes can lead to percolation, thereby forming an 

I infinite string network. This is remarkable as such an infinite string network is thought to arise only 

04 . when the entire universe undergoes phase transition. 

> : 

1-2 i Condensed matter physics experiments can provide very fruitful analogues of phenomena thought to have occurred 
in the early universe. An important example of this is provided by the realization that the theories of formation of 
0^ , exotic topological defects, such as cosmic strings in the early universe 1] can be tested in condensed matter systems [2|. 
' Many experimental studies have been carried out in liquid crystals [3] , superfluid helium 4} , and superconductors [5[ , 
^ ] by looking for analogies of cosmicphenomenon, as well as by focusing on universal properties for rigorous quantitative 
' tests of cosmic defect theories [a, in these systems. Here, we demonstrate the tremendous potential of this 
^f) \ correspondence by showing that a completely new insight into the phenomena in the early universe is obtained from 

■ a liquid crystal experiment, namely, the formation of large string loops due to convective stretching in a liquid crystal 
sample with a heated wire tip. We show that the same physics implies that local heating of plasma in the early universe 

^ I ' by evaporating primordial black holes [Si] can lead to formation of large cosmic string loops which can percolate to form 
\^ [ an infinite string network (even when ambient temperature is far below the cosmic string scale). This is remarkable 

■ as such an infinite string network is thought to arise only when the entire universe undergoes phase transition. 
Formation of cosmic strings is of great importance as they can provide a direct possible window to the physics 

^ • of ultra-high energy scales. Recent realization that superstring theories can also lead to existence of cosmic strings 
I I Pi has further added to the motivation of understanding various mechanisms of formation of cosmic strings. In the 

conventional picture, strings are produced when the universe goes through a symmetry breaking phase transition 
(-H This may not always be possible, for example, in models of inflation which have low reheat temperature, or in models 
' of TeV scale gravity [10|. Our results are therefore important as they provide a new possibility where an infinite 
. !^ ] string network can form without an over all phase transition and when the ambient temperature is much lower than 
^ • the cosmic string scale. 

H ] Our experimental setup consists of a droplet (size ~ 2-3 mm) of uniaxial nematic liquid crystal (NLC) 4'-Pentyl-4- 
- - 1 biphenyl-carbonitrile (98% pure, purchased from Aldrich Chem.) placed on a clean, untreated glass slide. For uniaxial 
nematic liquid crystals the order parameter is called the director which is a vector with opposite directions identified. 
This arises due to rod like nature of the molecules in this case (or, if the molecules are disk like, they are stacked to 
form a rod like structure) . In the low temperature phase (the nematic phase) , molecules tend to align locally, leading 
to non zero magnitude of the order parameter with the director representing the average orientation of the molecules 
in a local region. In the high temperature phase (the isotropic phase), which occurs above a critical value of the 
temperature, Tc, molecules are randomly oriented leading to vanishing magnitude of the order parameter. This leads 
to the order parameter space being RP^{= S"^ /Z2) which allows for string defects with winding 1/2. The director 
rotates by tt around this string defect. 

The isotropic-nematic (TN) transition temperature Tc for the NLC we have used is about 35.3 °C. The sample was 
kept at the ambient temperature less than Tc so it was in the nematic phase. A copper wire of about 100 micron 
thickness, with shape ' was connected to a resistor for heating (using voltage supply) such that the tip of the wire 
could be lowered into the sample drop from the top. The wire tip was in between the objective of the microscope 
and the sample. By increasing the current, the wire tip could be heated to a temperature much higher than Tc (upto 
about 50 - 60 " C). A Leica DMRM microscope with 20x objective, and a CCD camera was used for observations, 
at the Institute of Physics, Bhubaneswar. The string formation was recorded using a camcorder (for larger field of 
view) . 

When the current is increased with the wire tip inside the NLC drop then a small region around the tip heats up 
above Tc to the isotropic phase. When the current is low such that this isotropic region is small then small string 
loops form and shrink near the wire tip, as seen in Fig. la. As convection currents take the liquid (crystal) elements 
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away from the heated tip, the temperature of the Hquid element drops and falls below Tc at a certain distance from 
the tip. Resulting I-N phase transition can lead to string formation via the standard Kibble mechanism Another 
way strings can form here is because of turbulence generated by the heated tip due to temperature gradients in both 
vertical and horizontal directions. This can lead to formation of string defects as the order parameter (director) 
undergoes non-trivial rotations due to turbulent vorticity. 

Fig. Ib-d show large string loops further away from the wire tip as the wire tip is heated further. (Tip appears 
bigger in the picture when it is embedded deeper in the sample.) Larger heating increases string formation. At the 
same time we see that string loops, instead of shrinking, move away from the tip and stretch to large sizes. This 
is surprising, as one expects string loops in liquid crystal to always shrink due to string tension, (i.e. due to excess 
free energy). This is because, unlike superfluid vortices, there are no velocity fields associated with the winding of 
the liquid crystal string defects. What happens here is that strings are swept away from the wire tip by the radial 
convection currents which force strings to stretch and become larger. Eventually string loops shrinks, presumably 
when the fluid velocity subsides and is unable to compensate for the string tension (i.e., the tendency for string length 
to shrink due to excess free energy). Due to vertical temperature gradient, the convection currents also force string 
loops to return towards the wire tip. (Interestingly, confirmation for this pattern of convection currents is provided 
by observation of dirt particles embedded in the sample which are seen to follow the motion of string loops.) Often 
we see that string loops get trapped in cyclic motions, repeatedly swept away from the tip and becoming larger, then 
shrinking while returning towards the tip. We have checked that when the tip is kept just outside the sample (so 
there are no boundary effects due to the tip inside the sample) then string loops are still produced, though lesser in 
number (which, at least partly is due to lesser heating). 

There is a natural analog of this phenomenon in the early universe, where the Hawking evaporation of primordial 
black holes Q can lead to local heating of the ambient plasma 11, l3 | (se e, also, [l3|)i leading to restoration of various 
symmetries. An investigation of this local heating was carried out in [llj utilizing known results about the energy loss 
of quarks and gluons traversing a region of quark-gluon plasma [l3| . (We mention here that the results of ref. [lj| are 
certainly applicable for high energy partons (say several 100 GeV), it is not clear whether they can be extrapolated to 
the case of Hawking radiation with energies almost near the GUT scale. Still, the overall picture of pbh heating the 
loc al p lasma to very high temperatures may remain valid even for such large Hawking temperatures). It was shown 



ll[ that this energy loss leads to rapid heating of the plasma near the black hole. Further, resulting temperature 



gradients lead to large pressure gradients such that plasma develops radial flow. 

Possibility of formation of primordial black holes, and various observational constraints on them are well discussed in 
the literature, see refs. for early discussions. In first order phase transitions, particularly in first order inflationary 
scenarios, primordial black holes can be produced by collapsing regions of false vacuum ^15i] or due to inhomogeneities 
formed during bubble wall collisions |16||. Thev could also form by large amplitude density perturbations produced 
due to fiuctuations in the infiaton field |l7| , or by shrinking cosmic string loops fl^ . 

For primordial black holes produced from large amplitude density fluctuations produced by the inflation, one might 
ask whether CMBR data puts any constraints. However, as discussed in ref. [19|, one should note that CMBR data 
imposes restrictions on the power spectrum at large scales (of order 1 - 10^ Mpc), whereas, the primordial black holes 
form due to large amplitude and much smaller scale density fluctuations produced by the inflation. Hence they do 
not affect CMBR results at all but imposes constraints on the spectral index and consequently puts restrictions on 
certain models of inflation. 

The constraints on primordial black hole formation come from other sources like nucleosynthesis, gamma ray 
background etc. [l3| • However, all these phenomenon essentially imposes constraints on primordial black holes with a 
mass range which is much larger than the ones relevant here. We consider primordial black holes with masses at most 
few tonnes and such black holes evaporate away much above the QCD scale, without causing any conflict with current 
observations. Our main aim is to illustrate a completely novel possibility of forming inflnite string networks via back 
holes. Hence, we will simply assume the required masses and number densisties of black holes, without discussing any 
speciflc mechanism which could give rise to the formation of such black holes. 

Consider the case when the plasma heated by the black hole leads to local restoration of a gauge symmetry (with an 
energy scale and associated phase transition temperature of order 77) which allows for the existence of cosmic strings 
[l[. As the plasma flows out radially away from the black hole, its temperature will decrease and fall below Tc{— rj) at 
some distance where the symmetry will be broken. Near that region, again, string defects will form either via the 
Kibble mechanism, or due to turbulence. (For internal symmetries, spatial variations of the order parameter could 
lead to formation of non-trivial windings when turbulent motion of the plasma folds up extended spatial regions, 
essentially compactifying parts of spatial regions, i.e. plasma). If the friction forces dominate over the string tension 
so that strings are effectively frozen locally in the plasma then the string loops will be carried out by the radially 
expanding plasma. They will then stretch out to large sizes instead of shrinking due to tension, just as for the liquid 
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crystal experiment in Fig.l. It is important here to realize the difference between string tension for the liquid crystal 
example and for the relativistic cosmic string case. For liquid crystals, string tension represents increased free energy 
associated with the string defect configuration. String shrinks here as the free energy of the system decreases. For 
cosmic strings, it is the relativistic energy of the cosmic string configuration (arising from the effective potential) . 
String tension leads to shrinking of cosmic string when energy can be dissipated in other modes, such as gravitational 
waves, or dissipation in background plasma. The final effect is the same that the strings like to shrink. In our case, 
strings expand because plasma flow increases the energy (free energy for liquid crystals, and relativistic energy for 
cosmic strings). 

If there is a uniform density of similar primordial black holes, then each black hole will emit string loops which 
will be stretched to large sizes as they are carried away by convective flow, and may intersect each other. It is 
well established that string defects when crossing each other, intercommute [20| (unless they have ultra-relativistic 
velocities, which is not the case here). Intercommutation of string loops from different black holes will result in 
formation of larger string loops, leading to the possibility that these string loops may percolate and form an infinite 
string network. 

Let us discuss the conditions for percolation of strings. A black hole of mass Mbh evaporates by emitting Hawking 

radiation with an associated temperature Tth = st^Mbh '^^^^^ ~ ^-"^ ^ lO^^GeV is the Planck mass. We use 
natural units with H = c = 1. The rate of loss of mass by the evaporating black hole is given by 



dMbh _ aM^ 
dt ~ M? 



PI 



(1) 



bh 



Here, a accounts for the scattering of emitted particles by the curvature and depends on Tth- For values of Tbh 
relevant for our model, we take a to be ~ 3 x lO^'^ (see, [III El). The lifetime Tbh of the black hole obtained from 
Eq.(l) is Tbh — lO^Mp^Mg. where Mq is the initial mass of the black hole. Eq.(l) implies that very little energy is 
emitted until time of the order of Tbh which is when most of the energy gets emitted. The age of the Universe tu when 



its temperature is Tjj is [22| tu ~ 0.3g* ^^"^ MpiTjj^ , where (7*(— 100) is the number of relevant degrees of freedom 



|22l |. By equating Tbh = tu, we can get the mass Mq of the black hole such that its evaporation becomes effective 
when the temperature of the Universe is Tu, 



Mo = omM 



5/3^—2/3 
PI 



(2) 



For example, with Tu = 1 GeV we get Mq = 4 x 10"A/p;, Tbh = 10^ GeV and Uh = 5 x lO^''GeV-^ = 3 x IQ-"^ 
s. (We mention here that the mass of a primordial black hole can also increases due to accretion of background 
plasma particles. However, for relevant ranges of black hole masses here, all the dominant particle species are ultra 
relativistic, so only the geometric cross-section of black hole is relevant which is not very effective [23|. Some slow 
growth of the mass of the black hole could occur because of this, and the black hole masses we use here should be 
taken to be the final mass of the black hole when its evaporation becomes effective.) 



In ref. 11[, conditions for the equilibration of the Hawking radiation in the ambient plasma were discussed. Further, 



using the Euler equation for a relativistic fluid, it was shown in [11] that plasma velocity rapidly becomes relativistic 
due to large temperature gradients near the black hole. Assuming a steady state situation where the luminosity L{r) 
is independent of the distance r from black hole center, and is equal to —dMbh/dt (Eq.(l)), the temperature profile 
T{r) of the expanding plasma was obtained in 11| to be. 



r(r) = 



L 

(87r35,/45)r2i;p 



(3) 



This was valid for distances where bulk plasma flow dominates over diffusion of particles [IT] , and will be applicable 
for distances relevant here. We will take the plasma velocity Vp to be of order of sound velocity Vp ~ l/v^. We are 
neglecting shock formation here for simplicity, though shocks may be more favorable for string percolation. 

With Eq.(l) we see that it is the last stages of black hole evaporation which will lead to highest temperatures for the 
plasma in the nearby region, and will be most relevant for our model. Thus, we consider the situation when a black 
hole of initial mass Mq (Eq.(2)) has reduced by evaporation to a much smaller black hole with a mass M^ = Mq/x, 
whose life time is t^ and temperature is T^ > rj. Here x > 1, and is constrained by requiring that the thickness of the 
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heated plasma region (which is of the order of Tx for plasma expanding with sound velocity) should be greater than 
?7~^ (for self consistency of string formation at scale rf). We thus require. 

100M^-,^M > ,,-1 . (4) 

Temperature profile of the plasma around this remaining black hole is obtained from Eq.(3), and using Eq.(l) with 
Mbh = Mx = Even though steady state is not a good approximation for late stages of black hole evaporation, 
we will take this profile for our case. In our case it will not represent the entire temperature profile at a given time. 
Rather one can think of this as temperature of plasma shells which expand and cool. Assumption of constancy of 
luminosity may then be valid for these plasma shells. In this respect it is important to note that essentially the same 
profile is obtained if we take T = effectively at a distance of order of the Schwarzchild radius = xp- — along 

with the condition T"*r^ — constant (which essentially means conservation of energy-momentum of the expanding 
shells) . 

Temperature will drop to the value 77 at a distance which is obtained by setting r(r^) — rj for T(r) (obtained 
from Eq.(3)). String loops will be produced at this distance and will be dragged (and stretched) by expanding plasma. 
The plasma flow should stop at a distance of order ru where the temperature drops to the ambient value Tu. ru, 
therefore, sets an upper limit on the distance Rstretch up to which strings can be carried out by plasma flow, and 
is determined by setting T{ru) — Ty. Using Eqs.(2),(3), we get, ru ~ Q.{)AxM^[^T^'^^^ . To determine Rstretch we 
need to find when string motion is dominated by friction forces arising due to the scattering of plasma particles from 
the string. The dominant contribution to this comes from Aharanov-Bohm scattering (of particles with appropriate 
fractional charges) for gauge strings [l^l . The magnitude of the friction force per unit length on the cosmic string is 
then given by [2J], 

F/rtct~/3r3— ^ , (5) 
V 1 — w 

where v is the string velocity through the plasrna, T is the plasma temperature, and /3 is a numerical parameter 
related to the number of relevant particle species [24| . In specific grand unified theory models 25] it could be of order 
one, and we will assume that to be the case here. (See, also, ref. [2^ for early discussions of scattering of particles 
from strings.) We mention here that stretching of string due to plasma flow has been discussed by Chudnovsky and 
Vilcnkin in ref. [27| where they considered light superconducting cosmic strings getting stretched due to turbulent 
plasma flow in the galactic disc and even in stellar interiors. One could also consider formation and stretching of such 
strings around primordial black holes. This will require considerations of magnetic fields in the cosmic plasma etc. 

We consider a simple situation when a string encircling the black hole is stretched symmetrically by radial flow of 
the plasma. (Order of magnitudes of our estimates should remain same for other geometries of strings, e.g. string 
ends may lie on the boundary of the symmetry restored region around the black hole, as seen often in the liquid 
crystal case.) Such large strings should easily form during string formation on the surface of the sphere with radius 
r,, near the black hole (as also seen in Fig.l for the liquid crystal case). String tries to collapse due to its tension. For 



a string loop of radius R, one can estimate [27[ the tension force to be of order Ft^ ^, f± ~ H. 



tension ^ ■ 

At the formation stage, near r r^, the string will be at rest in the local (expanding) plasma frame with v — 
in Eq.(5), and hence Ff^ict = 0. String loop will then start collapsing under its tension (in the local rest frame 
of the plasma). This will lead to a non-zero and consequently a non-zero friction force via Eq.(5). String loop 
will expand in the rest frame of the black hole, due to drag of the plasma expanding with sound velocity, when 
Ffrictiv = I/V3) > Ftension- This Condition gives the upper limit for Rstretch as, 



Rstretch < 10 * X j^^pj^^ = Rmax ■ (6) 

As plasma flow stops beyond a distance of order ru, we conclude that plasma flow can stretch string loops to 
sizes of order Rstretch — 'min{Rmax, ru). We should also consider the effect of black hole gravity on the string loops. 
Gravitational force per unit length on the string due to the black hole of mass Mbh can be roughly estimated as '^^1''^' 
where R is the separation of the string segment from the black hole. (We are taking string as a simple gravitating 
system, as for a string loop. For special geometries, such as for a straighter string the gravitational force will be 
different, and will be typically smaller). We see that, per unit length, this gravitational force becomes much less 
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than the force fLension ^ ^ due to string tension when R > GMbh- Thus for any distances much larger than the 
Schwarzchild radius of the black hole, gravity of the black hole is subdominant compared to string tension forces. 
We are considering the situation when plasma drag forces completely dominate over the string tension forces, leading 
to stretching of string to large distances (where black hole gravity becomes even more negligible). Thus black hole 
gravity is not relevant for string dynamics in our case. 



TABLE I: Sample parameter values for percolation of strings 
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If black hole number density in the universe is such that inter-black hole separation dbh is smaller than RsLreich 
then loop intersections will be frequent. This will lead to intercommutation of strings with probability (almost) 1, 
implying percolation of strings, resulting in an infinite string network. By assuming that energy density of black holes 
contributes only a fraction / of the total energy density of the universe, we get (with ~ 100), 



dbh ^ . (7) 

Percolation of loops (forming an infinite string network) will happen when dbh < Rstretch ■ With Rstretch = I'u (when 
ru < Rmax), this gives the minimum allowed value of x as Xmin — 3/^^/^(Mp;/Tc/)^/^ (consequently, maximum value 
of Mx). We will give results with x = Xmin (so that dbh = Rstretch)- We find that for string percolation, with 
/ = 1,10^3 and 10""^, 77 must be below about 2 x lO", 10^2,2 x IQio, and 5 x 10^ GeV respectively, while 

Tu must be less than about 2 x 10^^, 10^,3 x 10^, and 5 x 10^ GeV respectively. (Note that should be large 
enough to create many strings of sizes ~ Rstretch- This is always true with our parameters.) When ru > Rmax then 
Rstretch = Rmax- V Can be larger in this case (e.g., up to 4 x 10^''' for / = 1). However, for very large 77, cannot be 
much bigger than Mpi. Use of Eq.(l) may be suspect for such small black holes. We will quote results for parameter 
values when ru < Rmax, and will not consider very large values of r]. 

In Table 1 we have given several sets of values for various parameters for which string percolation happens. We see 
that RsireLch » r^-, meaning that the strings are dragged (via friction forces) and stretched by the plasma flow up 
to very large distances. We also see that the percolation of strings is not a rare occurrence, but happens for rather 
generic set of parameter values. For example, even with almost negligible fraction of energy density in black holes, 
/ = 10"^, string percolation is achieved. 

We conclude by emphasizing that our results present first possible scenario where string loop generation by localized 
sources (heated wire tip for the liquid crystal experiment, and primordial black holes for the early universe) can lead 
to the formation of an extended string network without an overall phase transition, and even when the ambient 
temperature is much below the string scale. During scaling regime, there may not be any difference in the properties 
of string networks formed via the mechanism discussed here, and those arising from the conventional mechanism 
during a phase transition. These results illustrate the tremendous potential of correspondence between phenomenon 
observed in condensed matter systems and those which could have occurred in the early universe. 
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Figure Caption 



Fig.l: 

Production of string loops by heated wire tip embedded inside NLC droplet. Strings are carried out by convection, 
and are stretched to large sizes 



